Two-subband electron transport in nonideal quantum wells 
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Electron transport in nonideal quantum wells (QW) with large-scale variations of energy levels 
is studied when two subbands are occupied. Although the mean fluctuations of these two levels are 
screened by the in-plane redistribution of electrons, the energies of both levels remain nonuniform 
over the plane. The effect of random inhomogeneities on the classical transport is studied within 
the framework of a local response approach for weak disorder. Both short-range and small-angle 
scattering mechanisms are considered. Magnetotransport characteristics and the modulation of the 
effective conductivity by transverse voltage are evaluated for different kinds of confinement potentials 
(hard wall QW, parabolic QW, and stepped QW). 
PACS 73.20.Dx; 73.40.-c; 72.10.Fk 



> 1 I. INTRODUCTION 

co 
*n . 

The study of transport properties of semiconductor structures with many occupied electronic subbands has attracted 
CO • scientific interest for a long time, concerning especially with two-subband occupancy and screening effects in the scat- 
tering processes [0-0. Moreover, until recently these studies were made in inversion layers of Si-MOS (metal-oxide- 
semiconductor) structures in GaAs/AL r Gai_ ;r As selectively doped heterojunctions with two-subband occupied 
[^-^| , where large-scale fluctuations of layers have small effect on the position of the energy levels and transport prop- 
erties. Recently, new types of quantum wells (QW), e.g., modulation-doped QW @, wide parabolic QW §j||I), and 
stepped QW []l0| semiconductor structures with two occupied subbands, have been grown by tailoring the conduction- 
band edge of III-IV semiconductors. It was shown that these novel heterostuctures with double subband occupancy 
exhibit a much stronger effect of large-scale disorder on transport properties than the previously studied one-subband 
systems. For the new systems, the screening is effective only on the average of the external changes of the bottoms of 
the two occupied subbands, which are different from each other, and caused by smooth variations of QW parameters 
on the scale of the Bohr radius as illustrated in Fig. la. Only essential unscreened large-scale fluctuations of the both 
occupied QW levels are present in these systems. 

This paper addresses the effect of such non-screened variations of both occupied energy levels on the electronic 
spectrum and the classical magnetotransport by employing different confinement models: hard wall QW (Fig. lb), 
parabolic QW (Fig. lc), and stepped QW (Fig. Id). Short-range and small-angle elastic scattering mechanisms, within 
a local response approximation, are included in the linearized kinetic equation. We obtain effective magnetotransport 
coefficients, which are evaluated from the components of the local conductivity tensor averaged over large-scale 
inhomogeneities, up to second order in the disorder contributions to the conductivity {T^| , for zero temperature. In 
addition, we calculate the modulation of the effective conductivity by the transverse voltage. We take into account 
the difference between intra-subband momentum scattering frequencies and the inter-subband contributions. We 
investigate also the spatial variations of these frequencies. We find an essential broadening of the peak of intersubband 
transitions due to the large-scale inhomogeneities and we provide numerical estimates of magnetotransport coefficients 
for the considered models. 
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The organization of the paper is as follows. In Sec. II the energy spectrum, which takes into account non-screened 
variations of levels, and the general conductivity tensor within the linear response approach, are obtained. From 
the transition probability for a random potential model, we calculate, in Sec. Ill, the scattering frequencies and 
magnetotransport coefficients considering short-range and small-angle scattering mechanisms and after averaging the 
conductivity tensor over in-plane large-scale inhomogeneities. In Sec. IV, we discuss the assumptions, on which the 
work described here is explicitly dependent, and we present our concluding remarks. 



II. GENERAL EQUATIONS 

We discuss now our model for the energy levels fluctuations due to random inhomogeneities and we obtain the 
general expressions of the conductivity tensor based on the local response approach for nonideal QWs with two 
occupied subbands. 



A. Electron energy spectrum 

Electron energy levels in a QW with slowly varying width can be described in terms of r = {x, y} along with the 
in-plane kinetic energy given by the parabolic dispersion law e p = p 2 /2m with effective mass m. For instance, the 
influence of nonideal heterointerfaces in the hard-wall QW, formed by large band offsets, can be studied by assuming 
that the QW widths, d r , are randomly varied. Then, the few lowest energy levels s = 1, 2, ... are given by 

- -W^ 2 ~e 8 (l-2 S 4), (1) 



2m \ d 

where e s = s 2 (ttH/ 'd) 2 /2m, d is the mean width of the QW, and 5d r is its random fluctuation due to heterointerface 
roughness. The energy of an electron in the s-th subband is given as Epr — p 2 /2m + e sr + v sr , where v sr — 
J dz \ (fisz\ 2 V rz is the screening potential averaged over the s-th subband charge distribution \ip sz | 2 along the z direction. 
The potential V rz is determined by the Poisson's equation as 

\7 2 V rz = (n rz - (n TZ )), n rz = ^ n sr \tp sz \ 2 , (2) 

6 s=l,2 

where n sr — P2d(ef — £sr — v sr ) is the electron density in the s-th subband, p2D = m/'nh 2 is the 2D density of states, 
£f the Fermi energy, (n rz ) is the electron density averaged over the random inhomogeneities, and e is the dielectric 
constant, assumed independent of z. Hereafter, we consider only the case of two occupied subbands. 
The solution of Eq. (||) is given as 

V v = ^-J dz'e-^-^Sn^,, (3) 
where <5n qz is the 2D Fourier transform of the nonuniform part of the electron density 

Sn rz = ^2 \<Psz\ 2 5n sr = ~p 2D ^2\ ( p S z\ 2 (5£ sr + v sr )- (4) 

s=l,2 s 

Substituting Eq. (^) in Eq. (g), we are led to the system of two linear inhomogeneous equations for the Fourier 
transforms of screening potentials w sq for s — 1 and s = 2 subbands: 



27re< 
e<7 



^P2D J dz tf>z J dz'e'^- z '\ |yw | 2 (fe,,„ + tv q ). (5) 



Notice that inhomogeneity of Eq. (|5|) comes from terms oc fe S ' q . Hence the solutions are expressed in terms of the 
variations of the occupied levels, introduced by Eq. (Q). Assuming large-scale variations, £ c ^> d, it follows from Eq. 
(|J) that only the Fourier components with q < t^ 1 -C d^ 1 are essential. Then exp(— q\z — z'\) in Eq. (^) can be well 
approximated by the unity and we have 
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V S q = y^(<5£ s 'q + V s 'q), (6) 

s 

where is the Bohr radius. Here the right-hand side is independent of s, so i>i.2q = Vq. Further, for gas/8 <C 1, 
from Eq. (||), it follows that 



2 

8=1,2 

i.e., screening is relevant only on the average of the external variations fei r and 5e2v As a result, the electron 
spectrum energy can be written in the form 

E$ = ^+s s -{-iy5e r , (8) 

where Se r — (£2 — S\)5d r /d ss 3eiSd r /d is the non-screened part of potential. For further numerical estimates, we 
assume Gaussian correlations for the QW width fluctuations 

(6e r 6e r ,) = (6d r 6d r ,) = (fc) 2 w(\r - r'|). (9) 

Here 6e — 3e\6d/d, 8d is the typical height of the roughness interface and w(x) — exp[— (x/£) 2 } is the Gaussian 
function with lateral correlation length I. We can also write analogous expressions for the case of parabolic (Fig. 
lc) and stepped (Fig. Id) QWs but we do not give here the pertinent formulas considering further that 5s is a given 
quantity. 



B. Conductivity tensor 

Classical transport phenomena for electrons with spectrum, given by Eq. (@), are described by the distribution 
function fp(Epr > ) + AfpJ within the linear response theory. Here fp is the unperturbed Fermi distribution and the 
field-induced distribution is determined from the linearized kinetic equation of s-th subband given by 

(vV r + (-iyV6e r ■ A + f [v x H ] • A) A fg - 7 s (A/|pr) = -eE r • dfp jjf\ (10) 

where E r is the electric field, which is zero in the absence of a net current density j r , H is the magnetic field 
perpendicular to the QWs plane (v = p/m) and the collision integral for the elastic scattering case is written as 

/ s (A/|pr) = £ W^(A/^ - Afg). (11) 

s'p' 

VFpp, is the transition probability per second between states s,p and s',p'. The total current density j r is given by 

the standard expression j r = (2e/L 2 ) ^ sp vA/ P r' and satisfies the continuity equation V • j r = 0. 

The contributions from the convection and random fields in the left-hand side of Eq. (|l^) are estimated as vjl 
and Se/(p£) respectively. These contributions can be discarded when vf <C I and 6e e, where e = vp/2 is the 
characteristics electron energy, where we have used the relaxation time approximation to estimate the collision integral 
as r" 1 . Under these assumptions, Eq. ( [To|) assumes the form 

^[v x H] A A /W - J a (A/|pr) = -e(E r • v)/^ (12) 

where fp(E) = dfp(E)/dE. We look for a solution in the form 

A/pr' = ~ef' F (E$) {^(E r • v) + B s ^-([E r x H] ■ v)} , (13) 

where the coefficients A s and B s depend only on p 2 . Substituting Eq. ( fl3| ) into the kinetic equation (|i"2|), we obtain 
the linear algebraic system (see Ref. [0 for details) 
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B ° - E Kp' ( ~ C0S (P' p') ~ A *) = 1. 



s'p' 



^+E W pp' (^-cos(p,p')-S s J =0. (14) 

s'p' \ P / 

The solution of Eq. ( p^ ) for our two-level system is expressed through the transport relaxation frequency of s-th 
subband given by 

^^E^pp'^-^P'P')] (15) 
p' 

and the intersubband relaxation frequencies are written as 

E ll pp- = E M/ pp' cos (P'P')' ( 16 ) 
p' p' 

where s, s' = 1, 2. The solution of Eqs. ( O ) is 

Vl[u? + ^2 + J?1^12^2] - Vyi\p\zV<2, + T]i(uJ 2 + V 2 2 )] 



Ai = 



A 



o + ^2 + gjj + m^l2{v\ + ^2) 

^1 = "A ' ' ' 

with the cyclotron frequency w c = \e\H/mc, V\ = v l { + 1/12, 1^2 = + ui2, and A = [(^1^2 — v 2 2 ) 2 + co 2 {lo 2 + v\ + 
z/f + 2£f 2 )]- I n Eq. ( p7|) we have 771 = p'/p, where the momenta p' and p are connected through (p 2 — p' 2 )/(2m) = 
£2 — £\ — 25e r from energy conservation. The expressions for A2 and B2 follows from the above expressions for A\ 
and Bi, respectively, by changing 1 — ► 2, 2 — + 1 everywhere in Eq. ( ^T\) and taking 772 — p'/p where p' now satisfies 
(P 2 -p' 2 )/{2m)=e 1 -£2 + 25e r . 

Furthermore, for the zero-temperature, the coefficients A s and B s are calculated only at the Fermi level. Thus, we 
have to use p = pp 1 (p = pp 2 ) in Ai,Bx (A 2 ,B 2 ) while rji — Pf 2 /pf! and r) 2 = Pfi/pf 2 correspondingly. The induced 
current density, at zero temperature, is written as 

2 

j r = - V n« { A,E r + B s — [E r x H] I = a(r)E r , (18) 
to ' — ' l mc J 

s 

where the electron concentration for the s-th subband is given by 

4 s) =n s + (-iy P2D 6e T . (19) 



n s is the mean concentration of the s-th level and n = n\ + n 2 is the total average concentration. From Eq. (18), the 
components of the 2x2 conductivity tensor are 



p 2 . p 2 



fflI (r) = -V4 sIl ( s ), ^(r^-c^Vi^, (20) 

TO * — ' TO z — ' 

s s 

where cr^^r) = er y j,(r) and (Jj ;x (r) = —a xy (r). Note, that the conductivity dependence on r comes not only due to 

(s) 

the spatial variations of but also because pp 1 and pf 2 depends on r. 

III. TRANSPORT PROPERTIES 

In this section, we study the peculiarities of the steady-state magnetotransport caused by unscreened fluctuations 
of the QW's energy levels. In the model of random potentials, the transition probability can be written as 

W; s p , = ^Y, e ^-^/^) 2 }\{sV qz \s)\ 2 S(E p ^ s ,^ E p ^ r ), (21) 
9 
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where w exp(— Q 2 /q 2 ) is the Gaussian correlation function for the random potential with characteristic scale q 1 (note 
that go 3> Q = [(p' — p)/ft, q] is a 3D wave vector, and L is the normalization length. Summing over q in Eq. 

(|l|) leads to 

K ss , = I^ e -(^o) 2 | ( s '\e^\s) | 2 = ^ J dz j dz' Vs , zVszVs , z , Vsz ,e~^^/^ (22) 
and the transition probability transforms into 



W pp' = T72 exp{-[(p' - p)/hq } 2 }K ss ,S(E p ,^ :T - £ p , s , r ), (23) 



where the 5-function depends on r only for intersubband transitions, s ^ s'. Notice that Kyi = K 2 x- 

The functions K sa i are calculated below for the cases of short-range (q^ 1 -C d) and small- angle (hq <C Pf 

scattering. We consider three models of confinement potential: hard- wall QWs (see Sec. IIA), parabolic QWs 

with effective width £j_ = \J h/muj^, lo±_ is the characteristic frequency of the potential, and stepped QWs JTo| with 

quite different well widths d n <C d w . The coefficients K ss i are shown in Table I. 

Consider next the relaxation frequencies for these cases. For short-range scattering, straightforward calculations of 

Eqs. @, ©, and © give 

v[ r = (mw/h^Kn, v% = (mw/h 3 )K 22 , v 12 = (mw/h 3 )K 12 , (24) 
and i>i2 = 0. For small-angle scattering, we obtain 



5 « j?P§-K XXl and , 12 « f> 12 « - ^ ^W/^, (25) 

SV^PFia 2 V 7r « y / PF 1 PF 2 



where we assumed that (%o/pf s ) 2 "C 1 andpi? s is determined from £p^ = Ep. We point out that, in these conditions 
for small-angle scattering, we have V12 *C v 1 ^ due to the exponentially small factor. Indeed, for instance, in the case 
of the parabolic QW, V\ 2 Jv^y = \/pf 1 /pf 2 (pf 1 £±/K) 2 exp[— (jpp 1 — pf 2 ) 2 /ffq 2 ,]. For the maximum ratio pf 1 /pf 2 i we 
estimate vii/v^ — 4\/2 exp[— (p Fl — pf 2 ) 2 l^q®)- So, if exp[— (p Fl ~ Pf 2 ) 2 /h 2 qo] *C 10 _1 , we obtain v^/v*/ <C 1 for 
hard-wall and parabolic QWs. Similar condition holds for the stepped QW. For small-angle scattering, we neglect 



further very small contributions related to P12 and vi2- Then, as it is seen in the Table I and Eq. (25), the transport 
coefficients do not depend essentially on the specific QW for small-angle scattering. 

In order to analyze the influence of non- uniform contributions, we have to average the local conductivity tensor, 
Eqs. (po|), written in the form er M „(r) = {<j^ v ) + Xk(^°k)/ii/ exp(zk • r), where <5<Xk is the non-uniform contribution 
and {(Jfj,v) is the conductivity tensor for an ideal QW with two occupied subbands given by 

M= lj2_^\^ l* = x,v = x, ( } 

^ ' m ^ v I + u 2 { u c n = y, v = x 

Here we assume that intersubband frequencies are negligible. Following Ref. [Q, we obtain, up to the second order 
in the large-scale disorder, the effective conductivity tensor, from (j r ) = o e ^ (E r ), as 



eff - / \ _ ^«/3 ((^ cr -k)^fca fc /?(^ fJ k) ) 3 i .) 



where (...) denotes the average over the sample area, i.e., for distances r 3> £, and (SkC^O/w exp(ikr)) = 
p2[ . The average is performed using Fourier transform of the Gaussian correlation function, Eq. (0), (fe k fek') = 

5 k +k',o(7rfc 2 ^/ J L 2 )exph(fc4/2) 2 ] with k = Jkl + h*. 



A. Short-range scattering 

Before explicit calculation of the local magnetoconductivity tensor from the general Eq. (^0|), we present now 
standard expressions for magnetotransport coefficients of a QW with two occupied subbands. Using the conductivity 
tensor, given by Eq. (p6|), we obtain the magnetoresistance 
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Apj_ — ft — c 1 / — — TTT, (28) 

^1^2[(«1 + "2) + ("-1^2 + «2^l) J 

and the Hall coefficient 

R= 2_ ni ( V l+^) + n 2 ( V l+ul) 

ec [nf (z/f + w;J) + n|(f 2 + wjj) + 2nin 2 (^2^i + 

For strong (w c r ^> 1) and weak (w c r <C 1) magnetic fields, the Hall coefficient do not depend on H leading to 
R(oo) w l/ec(ni + "2) and = 0) ~ (nii^f + ri2i'i)/[ec(nxi^ + «2^i) 2 ]- 

For short-range scattering, all frequencies do not depend on Pf 12 and as a consequence, the coefficients A s and B s 
are independent of r. The non-uniform contribution to the conductivity tensor takes the form 

M v ; TO ^ z/2 + w 2 1 w c /x = y, V = X v 1 



After substituting Eq. ( |30[ ) into Eq. (|27]) and performing the averaging process, we obtain the components of the 
effective conductivity tensor as 

a eff- (a ) ^ f 1 1 1 gM + i^ ^ (31) 

^ xx; 2to 2 ( ffBX ) Uc 2 + »l "1 + z/f (c c 2 + ^ 2 ) 2 + w 2 (i/ 2 - Vl f ) ' 



e// _ / \ e4 Sn ( _L 1_\ / vi \ . , 

- to 2 ( ffM > Wc + *? - 2 + -ly U + -1 2 -c 2 + -i; ' ( j 

where <5rt = p2D$£ gives the concentration fluctuation due to large-scale inhomogeneities. 

Using Eqs. ( |3l| ) and (|32]), the effective magnetoresistance Ap e ff, and the Hall coefficient Reff, can be easily 
obtained. In Fig. 2, we plot Ap e // as a function of the magnetic field for n^/jii = 0.5. The solid and dot- 
dashed curves are 10 times enlarged and correspond to parabolic QWs where v^/vx = 5/6 and 5n/n\ = 0.2 and 
Sn/ni = respectively, The dotted and dashed curves correspond to stepped QWs for 5n/ni = 0.2 and 5n/ni = 0, 
respectively, with d n sw 10 A, d w sw 350 A, C/ w 0.1 eV, which leads to k ~ 10 6 cm -1 , £2 — £1 ~ 21.5 meV, and 
n-i +ri2 ~ 2.0 x 10 12 cm -2 . It follows that z^/^i = 1/2.1. Notice that Ap e ff = as well as Ap± = for the hard-wall 
QW because v^jv\ — 1 for short-range scattering. Then magnetoresistance is absent for hard- wall QWs, while a 
positive magnetoresistance appears for the parabolic QW, which becomes larger by a factor of 15 in the case of the 
stepped QW. 

In Fig. 3, we depict (R e ff — Rq)/Rq, where Rq = 1/ecn with n = (m + n 2 ) using the same curves and conditions 
as in Fig. 2. Similarly to Fig. 2, the solid and dot-dashed curves represent 10 x (R e ff ~ Ro)/Ro for the parabolic 
QW. In the case of the hard- wall QW, R e ff = Rq and R = i? because ^2/^1 = 1- Hence for hard- wall QW 
(R e ff — Ro)/Rq = 0, while an explicit dependence of H is observed for parabolic and stepped QWs. In addition, 
from Figs. 2,3 it follows that unscreened effect of the large-scale inhomogeneities of a QW parameters is manifested 
through stronger magnetoresistance at u> c /v-y < 1 and weaker at Lo c jv\ > 1. While it leads to smaller Hall coefficient 
for ijJ c jv\ < 0.5, for io c /v\ > 2 effect of such large-scale inhomogeneities on Hall coefficient becomes negligible. 
Notice that in experiment of Ref. || it was observed for magnetic fields, B z < 0.2 T, practically independent of H z 
magnetoresistance in GaAs-based parabolic QW with two occupied subbands: which leads to vi very close to v 2 that 
in turn can be valid only when scattering is short-range. 



B. Small-angle scattering 

For small-angle scattering, large-scale fluctuations of the electron density and scattering frequencies play an impor- 
tant role. As we have seen, scattering frequencies, given by Eq. (p5|), depend on fluctuations of the Fermi momenta. 
Further, we can neglect contributions coming from v\i and £12. Using Eqs. (|l7|), (20), (p5|), and (|27|), the components 
of the effective conductivity tensor, for small-angle scattering, can be written as 

-j— 2 

* e JJ/a [F 2 X (H) - F*(H)} , (33) 
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and 



rfjfo =5^- -^=F y (H)F x (H), 



(34) 



where 



2 £ T 



19 
16 



7 3gg/ig 
4 1 + wg/i^ 



(35) 



' yx 



= ((T yx ) /Wo = 2V2^2 



{u c /v s )(n s /n) 1/2 



8=1,2 



1 + LOI/V 2 



l + w c 2 /^ 2 



(36) 



and 



8=1,2 



5 _ 3gg/ij 
+ ^fV2 1 + c^f 



(37) 



- l^-D i + ^ I 4 - tt^m) ■ (38) 

Here z/ s = mwqQ/16Tfpp , where p F is independent of r and is determined by n s , ~5q — e 2 n/mvo, and a typical 

2/3 2/3 2/3 

scattering frequency vq is determined as v — [y x + v 2 )/2. Notice that za, is related to n s , while 
corresponds to n/2. In Figs. 4 and 5 we present our calculated effective conductivity tensor, given by Eqs. (|3^)-(^8|), 
magnetoresistance Ap e ff and the Hall coefficient, R e ff- 

In Fig. 4, we show Ap e ff for n-ijnx — 0.5 as a function of Lo c jv\. The solid and dashed curves represent Ap e ff 
for small-angle scattering for Sn/n — 1/8 and 6n/n — 0, respectively, and z^/fi = 8. For comparison, we plot 
the dotted curve for short-range scattering in the stepped QW as shown in Fig. 2. We observe in Fig. 4 that, 
for small-angle scattering, large-scale inhomogeneities strongly modifies the effective magnetoresistance in the case 
of two occupied subbands. Its behavior can be essentially different from the situation of short-range scattering. In 
this case, as it can be seen in Fig. 2, the magnetoresistance increases before reaching a plateau at uj c jv\ > 1. 
On the contrary, for small-angle scattering, this behavior is observed at oj c /vi > 8. As it is seen in Fig. 4, there 
is a much larger region of classical magnetic fields in which the magnetoresistance increases in case of small-angle 
scattering than for short-range scattering. Furthermore, in contrast with the short-range scattering case, we see that 
large-scale unscreened inhomogeneities lead to substantially smaller values of magnetoresistance at relatively weak 
magnetic fields, lo c /v\ < 3, while contributes to enhance the magnetoresistance at larger magnetic fields, io c jv\ > 4. 
Notice that according experimental results of Ref. p|, where all data are given for n.2/^1 ~ 1/2, the rather small 
magnetoresistance was attributed to small-angle scattering, which is confirmed by our calculated Ap e f f given by the 
solid curve in Fig. 4. 

Using the same parameters of Fig. 4, (i? e // — Ro)/Ro is depicted in Fig. 5. We now observe that large-scale un- 
screened inhomogeneities lead to a substantially weaker dependence of R e f f on H in agreement with the experimental 
results of Ref. |Tl[]. In conclusion, our theoretical results describe qualitatively well the experimental ones. jll| We 
point out that the results, in the case of small-angle scattering, are independent of QW confinement model. 



C. Transverse field effect on the conductivity 

Now we consider the modulation of the effective conductivity under transverse voltage which causes intersubband 
redistribution of the electron population. From Eq. (|3l|), for zero magnetic field, the effective conductivity, for 
short-range scattering, can be written as 

a eff W = 1 + — 1- 2— -3- — — - } , 39 

V nv ) 1 n 2 77 2 (l + AnAv nv) 2 \ 
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where a = e 2 nV jmvxv^, V = {y± + v%)/2, Av = {v 2 — v\)/2, An = (m — n 2 ) = Pid{^i — £1) with m = (n + An)/2 
and ri2 = (n — An)/2. For hard-wall QWs, because Av = 0, it follows that a e ^ is independent of An. In Fig. 6, we 
present Aa e ^ = a e ^ /a — 1 as a function of An/n. The dotted and dashed curves correspond to stepped QWs with 
Av/v ps —0.36, for the same parameters used in Fig. 2, and Sn/n =1/8 and Sn/n = 0, respectively. The thin-solid 
curve corresponds to a parabolic QW, with Avjv = —1/11, for Sn/n — 1/8 which practically coincides with that for 
Sn/n — 0. We observe that all curves, plotted in Fig. 6, for short-range scattering, have a negative slope. 
For small-angle scattering, the effective conductivity, given by Eq. (33), transforms, for H = 0, into 



r*ff l7f, 



/co 




An, 



5/2 



15Sn 
2n 2 



(1 



An, 



1/2 



(1 



An, 



1/2 



25<Sn 



[(i+^) 3/2 -(i-^f) 3/2 r 



2n 2 



{(l + 4r) 5/2 + (! - 4r) 5/2 + + 4r) 1/2 + C 1 - 4r) 1/2 ] } / 



(40) 



In Fig. 6, we show the results of Aa e ^ = a e ^ /ao — 1 as a function of An/n, calculated from Eq. (^) and represented 
by solid and dot-dot-dashed curves, which correspond to Sn/n = 1/8 and Sn/n — 0, respectively. We see that in 
contrast with the short-range scattering, now the curves have a positive slope. We do not discuss here the relation 
between An/n and the applied voltage. According to, for instance Refs. §, jy} and electro-optical measurements( 
for references see Ref. 



1 14 



the population redistribution can be essential already for rather small voltages. More 
detailed behavior can be achieved by performing self-consistent calculations for band structures with diagrams shown 
in Figs, l(b-d). 



IV. CONCLUDING REMARKS 



We have shown that non-screened variations of the intersubband energy modify essentially the transport properties 
of quantum wells when two subbands are occupied. Relative changes in the transport coefficients are determined 
by the degree of the random redistribution of the electron density and by the extent of the asymmetry of scattering 
processes that, for short-range scattering, are characterized by dimensionless parameters Sn/n and Av/V, respectively. 
Our model is founded on some assumptions, which we discuss below, and more detailed numerical calculations, which 
take into account peculiarities of the band structure and scattering mechanisms, are necessary for a full description 
of the electron transport. However, our results demonstrate that the kinetic coefficients change when we move from 
single to double subband occupancy in samples of similar quality. 

Along with magnetotransport characteristics (magnetoresistance, the Hall coefficient, etc.) and the field effect 
(modulation of effective conductivity by transverse voltage), other transport phenomena like Shubnikov-de Haas 
oscillations and cyclotron resonance are influenced by non-screened variations of subbands. For instance, essential 
broadening of the peak of intersubband transitions of a QW, due to unscreened large-scale inhomogeneities, |D| takes 
place also in the case of two occupied subbands. The effect of such inhomogeneities on intersubband optical transitions 
also requires special study. 

We discuss now the model assumptions. We have used the approximation of smooth inhomogeneities in several 
points: i) dependence of QW energy levels upon the in-plane coordinate r; ii) we have assumed that characteristic 
scale of smooth inhomogeneities i ^> as, and the screened potential, given by Eq. (0), results only on the average 
of variations of the energy levels; in) the condition that the transport scattering length vf <C t mainly determines 
the regime of applicability of the treatment; iv) the approximation of small amplitudes for smooth inhomogeneities, 
which leads to conditions Sn/nx,2 ^ lj reduces quantitatively the effects, but it allows to use theory applied to weakly 
large-scale inhomogeneous media fl^f . In addition, our calculations are based on simple QW spectra models instead 
of detailed self-consistent calculations of the band structure. However, we believe that the model assumptions, 
which are commonly used, do not change the overall behavior and numerical estimates of the magnetotransport 
coefficients. We hope our work can stimulate further experimental study of the transport properties of nonidcal QWs 
with multisubband occupancy in order to verify the effects of large-scale, unscreened inhomogeneities. 
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TABLE I. Coefficients K ss >, which appear in the transition probability given by Eq. (23). The columns labeled / 
and //stand for the short-range and small-angle scatterings, respectively, and for (6) hard- wall, (c) parabolic and (d) 
stepped QWs depicted in Fig. 1. hn — y / 2m{e2i — s w ) with e w = 2(ith/d w ) 2 /m, Kd n -C 1, nd w 3> 1, and q$d w <C 1 
and a w 4.6 x 1CT 3 . 





1(b) 


/f c ; 




//f 


11(c) 


77 W 




3/(2d) 


l/(V2vl±) 


k/2 


90/(20?) 


9o/(20?) 


9o/(2V^) 


K 22 


3/(2d) 


3/(4y/2ni±) 


3/(2^) 


90/(20?) 


9o/(20?) 


9o/(2V^) 


K l2 


1/d 


l/(2V2vr^) 


47r 2 /[d u ,(Kd u ,) 2 ] 


aq Q (q Q d) 2 


9oM±) 2 /(8Vi) 


327r 3 / 2 g 3/[ K 2 ( K d u; ) 3 ] 



FIGURE CAPTIONS 



FIG. 1. (a) Schematic views of the spatial variations of two occupied energy levels in the QW along the x direction 
with (solid curve) and without (dotted curve) screening. Band diagrams for (b) hard- wall, (c) parabolic and (cQ 
stepped QWs with non-ideal heterointerfaces. 



FIG. 2. Magnetoresistance Ap e // as a function of the magnetic field for short-range scattering. The solid and 
dash-dotted curves plot 10 x Ap e // in ca sc of parabolic QWs, Vijv\ = 5/6, with Sn/m — 0.2 and 8n/n\ = 0, 
respectively. The dotted and dashed curves correspond to stepped QWs, z^/^i = 1/2.1, with 5n/n\ — 0.2 and 
Sn/ni = 0, respectively. The ratio between average subband populations 712/711 = 1/2 and Sn is the typical value of 
the concentration fluctuations in both occupied subbands. 

FIG. 3. Hall coefficient R e ff as function of magnetic field for short-range scattering, where Z?o = 1/ecn, n = n 1 +n 2 - 
The curves represent the same cases as pertinent curves in Fig. 2. 



FIG. 4. Magnetoresistance Ap e // as function of magnetic field for small-angle scattering. The solid and dashed 
curves correspond to n<2,jn\ = 1/2 with Sn/n = 1/8 and Sn/n = 0, respectively. For comparison with the Ap e ff 
behavior for short-range scattering, the dotted curve is the same as in Fig. 2. 



FIG 5. Hall coefficient R e ff as function of magnetic field for small-angle scattering. The solid, dashed, and dotted 
curves correspond to pertinent curves in Fig. 4. 

FIG. 6. Conductivity Aa e ^ , for zero magnetic field, as a function of An = 711—712 showing the effect of intersubband 
population redistribution due to transverse electric field. The solid and dot-dot-dashed curves, with a positive slope, 
are for small- angle scattering and Sn/n — 1/8 and Sn/n — 0, respectively. The thin solid (parabolic QW with 
Sn/n = 0) dotted (stepped QW with Sn/n — 1/8), and dashed curves (stepped QW with Sn/n — 0) are for short- 
range scattering; other parameters are the same as for the stepped QW in Fig. 2. 
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